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Three-dimensional Voronoi analysis is used to quantify the clustering of inertial particles 
in homogeneous isotropic turbulence using data sets from numerics in the point particle 
limit and one experimental data set. We study the clustering behavior at different density 
ratios, particle response times (i.e. Stokes numbers St) and two Taylor-Reynolds numbers 
(Rex = 75 and 180). The Probability Density Functions (PDFs) of the Voronoi cell 
volumes of light and heavy particles show a different behavior from that of randomly 
distributed particles — i.e. fluid tracers — implying that clustering is present. The standard 
deviation of the PDF normalized by that of randomly distributed particles is used to 
quantify the clustering. The clustering for both light and heavy particles is stronger for 
the higher Re\. Light particles show maximum clustering for St around 1—2 for both 
Taylor-Reynolds numbers. The experimental dataset shows reasonable agreement with 
the numerical results. The results are consistent with previous investigations employing 
other approaches to quantify the clustering. We also present the joint PDFs of enstrophy 
and Voronoi volumes and their Lagrangian autocorrelations. The small Voronoi volumes 
of light particles correspond to regions of higher enstrophy than those of heavy particles, 
indicating that light particles cluster in higher vorticity regions. The Lagrangian temporal 
autocorrelation function of Voronoi volumes shows that the clustering of light particles 
lasts much longer than that of heavy or neutrally buoyant particles. Due to inertial 
effects arising from the density contrast with the surrounding liquid, light and heavy 
particles remain clustered for much longer times than the flow structures which cause 
the clustering. 



1. Introduction 

The distribution of particles transported by turbulent flows is a current resea r ch top ic 



with implication s in diverse fields, s uch a s process technology (IPratsinis et al 
cloud formation ( Bodenschatz et al. ( 2010h ). and plankton dynamics ( Schmitt fc Seuront 



cn top ic 
(Il996h) 



'. n most of the cases, the particles have a finite size and a different density than 
the carrier fluid, i.e. they have inertia. These inertial particles cannot totally follow 
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the fluid motion and distribute inhom ogeneously within the turbule nt flow, leading to 
clustering or preferential concentration (jToschi fc Bodenschatz ( 20091 )). The two relevant 
dimensionless parameters describing the dispersed inertial particles in the fluid are the 
density ratio j3 = 3pf/(pf + 2p p ), where pf and p p are the densities of the carrier fluid 
and particle, respectively, and the Stokes number, St= t p /t v , where t p = a 2 jZfiv is the 
particle relaxation time, t v is the typical timescale of the flow, which for a turbulent flow 
is the Kolmogorov time scale, a is the particle radius, and v is the kinematic viscosity of 
the fluid. 

In recent years, both numerical and experimental studies have quantified the cluster- 
ing of particles by employing different ap proaches like statistical a nalysis of sing l e-poin t 
measurement s (ICalzavarini et al. box- counting me t hod rtFessler et all J 1 994) ; 

Aliseda et all (|200j)), pai r correlati o n fun c tions (IChen et all (120061): ISaw et all t00$)). 



Kaplan - Yorke dimension ( Bee et all ( 20061 ) ; ICalzavarini et al I (12008^1. Minkowski func 
tionals (|Calzavarini et all ( 2008c )) and segregation indicators (|Calzavarini et al. ( 2008a[) ; 
IJzermans et ali i 20091 ) ). It is not possible to obtain global i nformation on bubble cluster 



ing from a single-point analysis ( Calzavarini et al. ( 20086f) ). Methods like box-counting 



and pair correlation functions, although useful, require the selection of an arbitrary length 
scale that affects the quantification of the clustering. The Kaplan- Yorke dimension, based 
on the calculation of the Lyapunov exponents, quantifies the contraction of a dynamical 
system by considering the separation rates of particle trajectories. Nevertheless, it does 
not provide global morphological information. Minkowski functionals, originally used to 
provide complete mor phological information of the large-scale distribution of galaxies 
(jKerscher et all (120011)). have been a pplied to study the cluste ring of particles in turbu- 
lent flows ( Calzavarini et al. ( 2008 ct )). ICalzavarini et al. ( 2008c ) found that light particles 
cluster in filamentary structures, whereas heavy particles have a wall-like topology around 
interconnected tunnels, and obviously no clustering was observed for neutrally buoyant 
tracers. In the above numerical simulations and experiments, the strongest clustering was 
found for particles with St«0(l). The problem with Minkowski-type analysis is that it is 
numerically expensive, and it does not provide information on the Lagrangian evolution 
of the clusters. 

An alternative mathematical tool that can be used to study clustering is the Vorono'i 
tessellation, which has been used in astr onomy as a tool to charact e rize c lustering 
of galaxies (jvan de Wevgaert fc Ickel (|l989l )). Recently, iMonchaux et all (|2010l ) apply a 
Vorono'i analysis to quantify the clustering of heavy particles in grid-generated turbu- 
lence. This Vorono'i approach does not require the selection of an arbitrary length scale 
for a fixed particle number, and it c an provide information on the Lagrangian statistics of 
clustering (jMonchaux et al. 2010l) . Monchaux et aZ| (l2010l) obtain two-dimensional par- 
ticle positions by imaging a turbulent flow in a wind tunnel seeded with droplets. The 
Vorono'i cells are defined based on the positions of the particles within the measure- 
ment domain. One can quantify the clustering by calculating the probability density 
function (PDF) of the normalized areas of the Vorono'i cells. The PDF will have a differ- 
ent shape for inertial particles when compared to the corresponding PDF of randomly 
distributed particles. The main difference is observed at the small and large values of 
normalized areas, where the PDF of heavy particles has a higher probability than for 
randomly distributed particles. There is a central region where there is no significant dif- 
ference between the PDFs of heavy particles and randomly distributed ones. The values 
of normalized areas at which the PDF deviates from the randomly distributed particles 
can be used as thresho lds to classify Vorono'i cells that belong either to clusters or voids. 
Monchaux et all (|2010l ) report a maximum preferential concentration for St around unity, 



in agreement with other methods that have been used to study clustering. 
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Table 1: Summary of the simulation and experimental parameters, where N is the size 
of the numerical domain, Rex is the Taylor-Reynolds number, r), t v are the Kolmogorov 
length and time scales, respectively, and N part i C ies is the number of particles in the 
simulations, and the time-averaged particle number in the measurement volume for the 
experiment. 



I N I Re x I n \ r v \ N par ucies \ St 



Simulation A 


128 


75 


0.0332 


0.1104 


1.0 x 10 3 


0.1 - 4.0 


Simulation B 


512 


180 


0.001 


0.0483 


6.4xl0 4 


0.1 - 4.1 


Experiment 




162 


288 /urn 


80 ms 


1.3xl0 3 


0.04 ± 0.02 



The objective of the present work is to extend the work of lMonchaux et al. I (|2010h to: (i) 
three-dimensions and (ii) a much larger range of density ratios (including light, heavy, and 
neutrally buoyant particles) and Stokes numbers, i.e. we quantify particle clustering by 
applying 3D Voronoi analysis both for numerical and experimental data sets of particles 
and bubbles. Moreover, we (iii) correlate the clustering behavior of different particles 
with local turbulent flow quantities and (iv) study the Lagrangian temporal evolution of 
the clusters. 



2. Experimental and Numerical Datasets and Voronoi analysis 

2.1. Datasets 

The numerical scheme for a dilute suspension (neglecting particle coll isions) of point 



partic l es in homogeneous and isot ropic turbulence is described as follows (jMaxev fc Rilev 
( 19831 ): ICaTzavarini et all (|2008ch ): 



where v = dx/dt is the particle velocity and u(x(t),t) the velocity field. The dimension- 
less numbers used to model the particle motion are the density difference between the 
particle and the fluid /3 and the Stokes number St. The values of j3 = 0, 1 and 3 correspond 
to very heavy particles, neutrally buoyant tracers, and bubbles in water, respectively. 
When St = 0, the particles perfectly follow the fluid flow behaving as fluid tracers. As 
summarized in table [TJ we explore a parameter space of /3 = 0, 1, and 3 and St ranging 
from 0.1 to 4.0 consisting of 24 values at Re\ — 75 with the spatial resolution of N = 128 3 . 
For Re x = 180 with N = 5 12 3 , we study 5 dif ferent values of St = 0.1 0.6, 1.6, 2.6, and 
4.1. (from iCFDdatabase [Fttp://cfd.cineca.it| ICalzavarini et alll2008c ). The simulation 



of the Navier-Stokes equation is based on a 2/3 de-aliased pseudo-spectral a lgorithm with 
2 nd order Adams-Bashforth time-stepping (for details see lBec et all (|2006[ )). Simulations 
have been performed in a cubic box of side L = 2ir with periodic boundary conditions. 
The forcing adopted acts only at the largest scale, it is implemented by keeping constant 
the kinetic energy content in of the smallest shell (\k\^ 1) in Fourier space. The intensity 
of the forcing is adjusted in such a way to have a turbulent dissipative scale (77) of about 
0.8 lattice grids in real space. Particle dynamics is evolved with time steps O(10) times 
smaller than the smallest Stokes time, leading to an accurate resolution of the particle 
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trajectories. Tri-linear interpolation is used to determine the value of the velocity field 
at the particle position. The numerical code was also validated by comparison against 
an independent code implementing different temp oral integration sche me, different parti- 
cles interpolation and different large scale forcing ([Toschi et al. The simulations 
extends over few 0(1) large-eddy-turnover times, this is enough for particles to reach a 
statistically steady distribution. In the present analysis, we fix the number of particles 
(Np ar ti C ies) for given Reynolds numbers: 1000 particles for the simulation with the do- 
main size of 128 3 at Re\ = 75, and 6.4xl0 4 particles for the simulation of 512 3 at Rex 
= 180. The number of particles normalized by the corresponding domain volume, i.e. the 
volume concentrations of the particles, for the two Re\ are identical. In one particular 
case of Re\ =75 and St = 0.6, the particle number is varied from 100 to 1 x 10 5 . 

We conduct experiments in the Twente Water Tunnel (TWT), an 8 m long vertical 
water tunnel designed for studying two-phase flows. By means of an active grid, nearly ho- 
mogeneous and isotropic turbulence with Re> up to 300 can be achieved. A measurement 
section with dimensions 2x0.45x0.45 m 3 with three glass walls provides optical access for 
the three-dimensional particle tracking velocimetry (PTV) system. Micro-bubbles with a 
mean radius of 170 ±60 /^m are generated by blowing pressurized air through a ceramic 
porous plate that is located in the upper part of the water tunnel. These micro-bubbles 
are advected downwards by the flow passing through the measurement section. In our 
3D-PTV micro-bubble experiments, we use a 4-camera system to get micro-bubble po- 
sitions in the active-grid-generated turbulence in the TWT. The experimental data are 
collected for a duration of 6 seconds (three times the large eddy turnover time) at an 
acquisition rate of 1,000 fps. For the experimental data, Re\ — 162, j3 — 3 and St = 
0.04 ± 0.02, and the time-averaged number of particles inside the measurement volume 



of 70mm x 7 0mm x 70mm is 1.3 x 10 3 (for further details, see iMartinez-Mercado et al 
(|20lll201ll) V 



2.2. Voronoi analysis 

The Voronoi diagram is a spatial tesselation where each Voronoi cell is defined at t he 
particle location based on the distance to the neighboring particles (jOkabe et al. I (|2000l) 1. 



Every point in a Voronoi cell is closest to the particle position compared to the neigh- 
boring particles, the exceptions are the vertices, borderlines and facets (see figure [T]). 
Therefore, in regions where particles cluster, the volume of the Voronoi cells is smaller 
compared to that of the cells in neighboring regions. Hence, the volume of the Voronoi 
cells is inversely proportional to the local particle concentration. The PDF of the Voronoi 
volumes normalized by the m ean volume for random ly distributed particles can be well 
described by a r-distribution (jFerenc fc Nedal (|2007h l (see figure [2]). In the 3D case, the 



r-distribution has the following prefactor and exponent: 

31 25 

f(x) = —x'e X p(-5x). (2-2) 

Here x is the Voronoi volume normalized by the mean volume. Particles which are not 
randomly distributed will have a PDF that deviates from this r-distribution, indicating 
preferential concentration. The Voronoi cells of particles located near the edges of the 
domain are ill-defined, i.e. they either do not close or close at points outside of the 
domain. These cells at the border of the domain are not considered for the analysis. 



Three-dimensional Lagrangian Voronoi analysis for clustering 



5 




Figure 1: An example of a 3D Voronoi tesselation. The dots represent the particle position 
and lines represent the borders of the Voronoi cells. 



3. Results 

First, we present results on the effect of the density ratio (/3) on the clustering, followed 
by the effect of the Stokes number (St) and the number of particles (N partic i es ) . Then, we 
show how the volume of Voronoi cells (V) and enstrophy are related. Finally, we present 
results on the Lagrangian autocorrelations of Voronoi volumes and enstrophy. 



3.1. Density effect 

Here we study the clustering behavior of particles of different f3 at a fixed St for two 
different Re\. Figure shows the PDFs of the Voronoi volumes (V) normalized by their 
averaged volume (V), V/V, for heavy, neutrally buoyant, and light particles of St = 0.6 
at Re x =75 (Fig. [2] (a)), and 180 (Fig. [2] (b)). It clearly shows that the trends in the 
probability density functions are similar for both Re\. The PDF of neutrally buoyant 
particles follows the T-distribution eq. (|2.2j) quite well, reflecting that neutrally buoyant 
particles do not have any preferential concentration. In contrast, the PDFs of light and 
heavy particles clearly show a different behavior compared to the randomly distributed 
particles. We observe that the probability of finding either small or large Voronoi volumes 
is higher for both light and heavy particles. The two regions of small and large volumes 
can be used to identify clusters and voids. The strongest clustering is observed for light 
particles, as the probability of finding small Voronoi volumes is the highest. Owing to the 
density difference , light particles accumulate i n vor te x filaments du e to ce ntrifugal forces 
(jMazzitelli et all (j2003l) : iMazzitelli fc Lohsel (12004) : iBiferale et all (|2010l )). while heavy 



particles concentrate in regions of intense strain ( Bee et al. ( 20061) ). Here, although the 



heavy particles show clustering, it is less compared to light particle s. These results are 
consistent with the Minkowski analysis by ICalzavarini et all (|2008a ) . 



3.2. Stokes number effect 

In this section, we study the effect of St on the clustering behavior for the three types of 
particles. We study the clustering behavior of the particles by examining the deviations 
of their Voronoi volume PDFs from the T— distribution. 

Figure [3] shows PDFs of light (/3=3) and heavy (/3=0) particles for different St at Re\ 
of 75 and 180. Firstly, we discuss the clustering of light particles as shown in Fig. [3] for 
Re\ of (a) 75 and (b) 180. Both Reynolds numbers give a similar trend with increasing 
St. When St increases, the probability of finding clusters and voids increases up to a value 
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Figure 2: (Color online) The normalized Vorono'i volume PDFs for heavy (squares), neu- 
trally buoyant (circles), and light particles (triangles) at St = 0.6 from DNS at (a) 
Re\ — 75, and (b) Re\ = 18 0. The thick lin e show s the T-distribution eq. (|2.2j) for ran- 
domly distributed particles ( Ferenc fc Nedal (|2007h ). the PDF of the neutrally buoyant 
particles agrees well with the randomly distributed particles (+). Both heavy and light 
particles show clustering, however light particles show the maximum clustering. 





Figure 3: (Color online) The normalized Voronoi volume PDFs for different St ranging 
from 0.1 to 4 in the numerics for (a) light particles /3=3 at Re\ = 75, (b) light particles 
(3=3 at Re\ = 180, (c) heavy particles /3=0 at Re\ = 75, and (d) heavy particles /3=0 
at Re\ = 180 . The stars in (b) correspond to the experimental result with St = 0.04 ± 
0.02 at Re x = 162. 
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Figure 4: (Color online) The comparison between log-normal and T— distribution of 
Eq. 13.11 fitting for the PDF of the 3D Voronoi volumes. Open symbols represent heavy 
(squares), neutrally buoyant (circles), and light particles (triangles) at St = 0.6 from 
DNS at Re\ — 75; the lines represent (a) log-normal, and (b) T— distribution. 



of St = 1.6, after which the dependence becomes weaker for both Re\. We note that the 
experimental result, shown with stars in Fig. [3](b), for micro-bubbles with St = 0.04 ± 
0.02 agrees reasonably well with the trend of the numerical data for light particles. In any 
case, for these small Stokes numbers, the PDF of the Voronoi volumes is still qualitatively 
similar to that of tracers. Another important feature of the light particle PDF is that 
the highest probability occurs at the smallest volume and de creases monotonically wit h 
increasing volume for St in the range of 0.6 to 4. As studied bv lCalzavarini e£~al ( 2008c ). 



bubbles in this range of St tend to get trapped in vortex filaments, leaving void regions. 
Thus, most of the bubbles are concentrated in small regions and there are few bubbles 
outside these small regions. 

In general, the clustering of heavy particles is weaker as compared to that of light 
particles. For heavy particles, as shown in Fig. [3] (c,d), as St increases, the probability 
of finding clusters and voids increases up to a value of St = 1.6, then the St dependence 
c hanges for different Re\ and is discussed below. 

Monchaux et al. I (|2010l) found that the Voronoi area statistics of heavy particles can 



be well fitted by a log-normal distribution. For comparison, figure 2] (a) shows that the 
log-normal fitting for the PDFs of the present 3D Voronoi volumes for the three different 
type of particles. It is clear that the Voronoi volume statistics cannot be characterized 
well by the log-normal function, even not for the neutrally buoyant particles. However, 
the PDFs for all type of the particles can be fitted very well by the T— distribution (see 
figure H] (b)) with only one fitting parameter a: 

/Or) = 2 1 - x^ 1 exp-^ (3.1) 

^r(^) 

where, a is the standard deviation of the Voronoi volumes. Hence a provides a proper 
statistical quantification of Voronoi volumes. 

In order to quantify the clustering using a single number, we use the standard deviation 
<t of the normalized Voronoi volume distributions. In figure[5](a), we plot a normalized by 
the standard deviation of the Voronoi volumes for randomly distributed particles or- The 
magnitude of the indicator a/or distinguishes the behavior of light, neutrally buoyant, 
and heavy particles. A higher value of the indicator reflects stronger clustering for a given 
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Figure 5: (Color online) (a) Normalized standard deviation a/ar (indicator) of the 
Voronoi volume distributions versus St for the two different Re\ from DNS data. The 
symbols correspond to heavy (squares), neutrally buoyant (circles), and light particles 
(triangles). Open and filled symbols represent data at Re\ — 75 and Re\ = 180, respec- 
tively. The value of the indicator for neutrally buoyant particles remains constant at 1, 
i.e. clustering is not observed, whereas light particles show the most clustering with a 
peak at St « 1.5 for both Re\. The experimental result of micro-bubbles is plotted with 
the star, (b) An enlarged plot showing only the results for light particles. 



i?e a- For neutrally buoyant particles there is no observed clustering, hence the indicator 
value is constant at 1. Heavy particles show clustering and the indicator value saturates 
at St ?s 1-2 at Re\ =75. However, the indicator value continuously increases with St at 
the higher Reynolds number of Re\ =180, and the absolute value of the indicator a/ar is 
larger for higher Re\ . This indicates that the clustering of heavy particles is stronger at 
higher Re\ for a given St. Fig. [5] (a) shows that the absolute value of the indicator a/ar 
for light particles is also larger for higher Re\, revealing a stronger clustering for light 
particles at higher Re\. The reason for the Reynolds number effect could be because of 
the changing range of length scales of the vortex filaments which affect the clustering. 
At higher Re\, there is a wider range of clustering length scales resulting in a Voronoi 
volume distribution with a higher value of standard deviation. The curves corresponding 
to light particles show the strongest clustering, with a peak at St ~ 1 — 2 for both Re\ = 
75 and 1 80. This clustering result has a consistent trend with that of the Kaplan- Yorke 
analysis ( Calzavarini et al\ ( 2008(J) ). 



We also add the data point for the standard deviation of the experimental Voronoi 
volume PDF as shown in figure [5J Although the mean value of the indicator a/ar for 
the experimental data is higher than those from the numerical simulations of light point 
particles, there is a good agreement with numerical trend within the experimental error- 
bar. More experimental data at larger Stokes numbers, i.e., larger bubbles, have to be 
taken to come to a final conclusion on this issue. 

3.3. Effects of the number of particles 

In principle, one can expect different behaviors depending on the number density of 
particles. In the simulation dataset A, there are 10 5 particles available for one special 
case of Re\ = 75 and St = 0.6. Using this snapshot, we study the effects of the particle 
number on the value of the clustering indicator a/ar- We subsample data from this 
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Figure 6: (Color online) Normalized standard deviation of the Voronoi volume distribu- 
tions as a function of number of particles taken from a snapshot case at Re\ = 75. The 
dashed line shows the number of particles we used in the present work at Re\ — 75 and 
St = 0.6 (simulation A). 



snapshot by selecting the required number particles and computing the Voronoi statistics. 
This subsampling procedure is randomized and then carried out at least 100 times for 
each case of particle number. Figure [5] shows the effect of varying the number of particles 
on the clustering indicator o~/o~r and the error bars represent the standard deviation 
of all the subsamples of a given number of particles. In the present data set, the mean 
distances of particles are: 34.477? for N part icies—^0 2 , I677 for N part icies—^0 3 , 3.44?y for 
N P articies=W 5 , which are all above In. Hence, we are always studying situations where 
the mean particle distances are in the inertial range. 

As shown in figure El for light and heavy particles the value of the indicator increases 
as the number of the particles is increased. The evolution of the value of the indicator 
a /or is steeper with increasing number of particles, and there seems to be no plateau 
region where the indicator value saturates. We do not understand the exact reason for 
this particle number d ependence. One possible reason could be that the clusters have a 



complicated structure ([Calzavarini et aLll2008d) . However, for a given number of parti- 



cles, the indicator does show a consistent trend: a stronger clustering for light particles, 
weaker clustering for heavy particles, and no clustering for neutrally buoyant particles. 
Moreover, the error of the indicator calculated at N part i C ie = 1000 is less than 4 %. 
Therefore, at a fixed number of particles, the clustering indicator a far of the Voronoi 
volume is robust. In the analysis that follows, we use the data of N part i C ies = 1000 for 
the simulation with the domain size of N = 128 3 at Rex = 75 (simulation A). 

3.4. Relation between the volume of the Voronoi cell and enstrophy 

We relate the Voronoi volumes for the three different type of particles with turbulent flow 
quantities. A natural property for t his co mparison would be the enstrophy £1 = lj 2 /2 
(where u is vorticity). IBenzi et~al. ( 20091 ) have shown that different types of particles 
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react sensitively to the local enstrophy at the particle position, reflecting their tendency 
to stay in regions with different vorticity contents. We thus calculate the joint PDF of 
Voronoi volumes and enstrophy for three types of particles at a fixed St = 0.6 for Re\ = 
75. For comparison, we also calculate the joint PDF for the case of neutrally buoyant 
particles with the smallest St available in the simulations (St = 0.1 and — 1). The 
statistical behavior of these particles is expected to be close to that of ideal fluid tracers 
(St = and f3 = 1). From now on, we refer to this case as the fluid tracer case. The 
Voronoi volume and the enstrophy are normalized by the mean values (Vtr and Qtr) of the 
fluid tracers. Figure [7] shows the joint PDFs of the normalized Voronoi volume (V /Vtr) 
and the normalized enstrophy (f2/f2 tr ) for the different type of particles. The joint PDF 
for neutrally buoyant particles of St = 0.6, shown in figure [7] (c), is very similar to that 
of fluid tracers shown in figure [7] (a). We observe a clear difference in the joint PDF 
for heavy and light particles, as shown in figure [7] (b, d). The coordinates corresponding 
to the peak of the joint PDF ((n/n tr )f^ F , (V p /Vi r )$?g F ) is indicated by the crosses 
in the figure for each case. Compared to the tracer case, a slightly lower (V p /Vtr)jpDF 
and a lower (Q/Q t r)JpDF f° r heavy particles indicates more clustering at low enstrophy 
regions. The maximum value of the joint PDFs for the light particles is located at the 
region with a much higher enstrophy and a smaller Voronoi volume. This shows that the 
light particles shows strong clustering at high enstrophy regions. 

The St dependence on the peak coordinates of the joint PDF ((0/O tr )^ F , (V p /Vtr)™PD F ) 
is plotted in figure 13 As shown in figure [8] (a) , the value of (V p /Vtr)JpDF f° r neutrally 
buoyant particles is nearly same with that of tracers at St from 0.1 to 4. The value 
of (V p /Vtr)Jp£> F f° r heavy particles is slightly smaller than unity for all St, indicating 
clustering. Figure [8] (a) also shows that the clustering for light particles is stronger as 
evidenced by the much smaller (V p /Vtr)Jpi>F compared to those of neutrally buoyant 
and heavy particles at all St. The minimum value of (V p /Vtr)Jp£>F indicating strongest 
clustering for the light particles is located at St = 1 — 2, which is in excellent agreement 
with the results obtained using the indicator er/or (figure [5]). The corresponding enstro- 
phy at the peak ((^/^tr)^pEF) °f the joint PDF versus St for the different particles is 
shown in figure [5] (b) . The value of (fl/n tr )^?g F for the heavy particles is smaller than 
unity, and it is much larger than unity for the light particles. This reflects the clustering 
of light particles in flow regions with very high enstrophy, whereas heavy particles cluster 
in low enstrophy regions for all St in the present study. 

3.5. Voronoi' Lagrangian autocorrelation 

Finally, we conduct a Lagrangian analysis on the Voronoi volumes. For each type of 
particle we calculate the Lagrangian autocorrelation of its associated Voronoi volume. 
Figure |9] (a) shows a typical temporal evolution of Voronoi volumes for the three types 
of particles at St = 0.6 and Re\ = 75. To compare the behavior of the three different 
particles, we choose particles with similar Voronoi volume at the starting time and trace 
their time evolution. While the Voronoi volumes of heavy and neutrally buoyant particles 
change frequently in time, it is clearly seen that light particles tend to have small values 
for longer times. This suggests that light particles are trapped in clustered regions for a 
long time and are suddenly ejected, as seen in figure [9] (a) around t/t^ « 95. 

Figure E] (b) shows the autocorrelation function Cy(r) for heavy, neutrally buoyant, 
and light particles at a fixed St = 0.6 and Re\ = 75. We define the decorrelation time ry 
as the time when the autocorrelation function has decreased to 1/2, i.e., Cy(ry) = 1/2. 
As shown in figured (b), the decorrelation time for light particles is around ry ~ 7t v , 
whereas for heavy and neutrally buoyant particles decorrelation already occurs around 
4t v . Thus the clustering of light particles lasts for a longer time as compared to heavy 
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Figure 7: (Color online) Joint PDFs of normalized Voronoi volumes and enstrophy for 
tracers and particles at St = 0.6 for Re\ = 75: (a) fluid tracers, (b) heavy particles, (c) 
neutrally buoyant particles, and (d) light particles. The cross indicates the location of 
the maximum probability (peak) for each case. 




Figure 8: (Color online) The coordinates of the peak of the Joint PDFs of normalized 
Voronoi volumes and enstrophy as a function of St: (a) V p /Vtr versus St, (b) f2 p /f2 tr 
versus St. 
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and neutrally buoyant particles. As shown bv ICalzavarini et al. ( 2008ch . light particles 
accumulate in filamentary structures and heavy ones tend to cluster outside these struc- 
tures to form wall-like interconnected tunnels. These differences in the morphology of 
the clustered particles could be a possible reason for the light particles being clustered 
for a longer time as compared to heavy particles. 

We also compare the autocorrelation time scale of the Voronoi volumes to that of 
the enstrophy shown in figure [H^c) for the same St and Re\. First, as expected, for 
neutrally buoyant particles, the Lagrangian decorrelation time for the Voronoi volumes 
is comparable to that of the enstrophy (to), i.e. to ~ Ty ~ &t v , because the neutrally 
buoyant particles do not cluster. However, remarkably, for light particles the decorrelation 
time of the Voronoi volumes is much larger, Ty ~ 7t v , i.e., more than twice as large as 
the autocorrelation time scale to ~ 3t v of the enstrophy itself. For heavy particles, the 
Lagrangian decorrelation time of the Voronoi volumes is around Ty ~ 4r^, which is also 
about two times that of enstrophy to ~ 2t v . 

We also study the St dependence of the decorrelation time scales of Voronoi volume 
(ry) and enstrophy (to) at Re\ = 75 for heavy, neutrally buoyant, and light particles as 
shown in figure [TU] (a). We observe that tv for light particles is always larger than heavy 
and neutrally buoyant particles in the St range 0.1 to 4, with a peak around St unity. 
This suggests that the light particles cluster for a longer time in the range of studied St. 

It is well known that flow regions of high enstrophy trap bubbles and regions with 
intense strain accumulate heavy particles. Figure ITOTa) also shows that Ty for both light 
and heavy particles is much larger than their decorrelation time of enstrophy to for all 
St from 0.1 to 4. This is more clearly seen in figure ITOT b) where the ratio tv/tq for both 
light and heavy particles is greater than unity for all St, while this ratio is always close to 
unity for the neutrally buoyant particles. Remarkably, this means that the life-times of the 
clustered bubbles and heavy particles are much longer than the life-time of the trapping 
flow structures themselves. The interpretation is that clustered particles are constrained 
in different regions of the flow and due to their inertia need time to reorganize themselves 
in the flow after sudden changes in flow conditions. However, neutrally buoyant particles 
do not have this constraint and are distributed more evenly at any given time in the flow. 
Figure [TUT b) shows that the ratio tv/tq for light particles has a weakly decreasing trend 
at St larger than unity. The ratio Ty /to for heavy particles monotonically increases with 
increasing St, and it is larger compared to light particles for St > 0.5. 



4. Conclusion 

We use three-dimensional Voronoi analysis to study particle clustering in homogenous 
isotropic turbulence with both numerical data in the point particle limit and one exper- 
imental data set. The analysis is applied to inertial particles (light, neutrally buoyant, 
and heavy) of different density ratios (3, St ranging from 0.1 to 4 and two different Taylor- 
Reynolds numbers (Re\ = 75 and 180). In the entire range of parameters covered, the 
Voronoi volume PDFs of neutrally buoyant particles agree well with the r-distribution 
for randomly distributed particles. At a fixed value of St, the PDFs of Voronoi volumes of 
light and heavy particles show higher probability to have small and large Voronoi volumes 
than randomly distributed particles, reflecting the clustering behavior. The standard de- 
viation of normalized Voronoi volumes ct/ot is used as an indicator to quantify the clus- 
tering. Heavy particles show some clustering, and light particles have a much stronger 
clustering. Both heavy and light particles show a stronger clustering for a higher Re\. 
The maximum clustering for light particles is around St w 1-2 for both Taylor-Reynolds 
numbers, and this maximum clustering range has a consistent trend with that of the 
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Figure 9: (Color online) Lagrangian Voronoi analysis for heavy (squares), neutrally buoy- 
ant (circles), and light particles (triangles) at St = 0.6 and Re\ = 75. (a) Temporal evo- 
lution of Voronoi volumes, (b) Temporal autocorrelation functions of Voronoi volumes, 
(c) Temporal autocorrelation functions of enstrophy. 

Kaplan- Yorke analysis. We check the effect of number of particles on the value of the 
indicator and find that the clustering trend is robust for a given number of particles. 

For one (small) Stokes number St = 0.04 ± 0.02 we have also extracted the 3D Voronoi 
volume PDF from experimental data. Though the PDF fits into the general trend - at 
these small Stokes numbers the PDF nearly follows a T-distribution - a quantitative 
analysis shows that the experimental PDF of 3D Voronoi volumes is slightly broader 
than what is obtained from point-particle simulations. More experiments with larger 
Stokes numbers will have to be done to judge whether this is a limitation of the point- 
particle approach, a consequence of the neglectance of two-way and four- way coupling in 
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Figure 10: (Color online) (a) Decorrelation time of Voronoi volume (ry) and enstrophy 
(rn) as a function of St at Re\ — 75 for heavy (squares), neutrally buoyant (circles), 
and light particles (triangles). Open and filled symbols represent decorrelation time of 
Voronoi volume and enstrophy, respectively, (b) The ratio of decorrelation times (ry / 
rn) as a function of St. 

the numerics, or whether the experimental data are not precise enough. From our point 
of view, the Voronoi analysis is an excellent means to quantitatively compare clustering 
effects of particles in experimental and numerical data sets. 

Finally, we show that the Voronoi analysis can be connected to local flow properties 
like enstrophy. By comparing the joint PDFs of enstrophy and Voronoi volumes and their 
Lagrangian autocorrelation, the clustering behavior of heavy, neutrally buoyant, and light 
particles can further be distinguished. It is found that the light particles strongly cluster 
in flow regions with very high enstrophy, whereas heavy particles weakly cluster in low 
enstrophy regions for all St in the present study. From the Lagrangian autocorrelation 
of Voronoi volumes we conclude that the clustering of light particles lasts much longer 
than that of heavy or neutrally buoyant particles. And due to inertial effects owing to 
the density difference from the carrying fluid, light and heavy particles remain clustered 
for a much longer time than the flow structures themselves. 
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